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Abstract 

By using a variant of the quantum inverse scattering method, commutation relations 
between all elements of the quantum monodromy matrix of bosonic Massive Thirring 

(BMT) model are obtained. Using those relations, the quantum integrability of BMT 
0\ '. 

model is established and the 5*-matrix of two-body scattering between the corresponding 
^O 

quasi particles has been obtained. It is observed that for some special values of the 

coupling constant, there exists an upper bound on the number of quasi-particles that can 
form a quantum-soliton state of BMT model. We also calculate the binding energy for a 
*^h. iV-soliton state of quantum BMT model. 
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1 Introduction 

Quantum integrable field models in 1+1 dimensions are objects of interest due to their 
close connections with different areas of physics as well as mathematics [1-10]. These inte- 
grable theories have played an important role in understanding the basic nonperturbative 
aspects of physical theories relevant in the realistic 3+1 dimensional models. Through 
quantum inverse scattering method (QISM) one can establish the integrability property 
of these models and obtain the spectrum as well as different correlation functions of the 
corresponding models [4]. 

Massive Thirring model in 1+1 dimensions has been widely studied as a toy counter- 
part to low energy QCD, since it does not include many of the complications arising in 
3+1 dimensions. The study of a nonlocal massless Thirring model is relevant, not only 
from a purely field theoretical point of view but also because of its connection with the 
physics of strongly correlated systems in one spatial dimension. This model describes an 
ensemble of non-relativistic particles coupled through a 2-body forward-scattering poten- 
tial and displays Luttinger-liquid behaviour [11] that can play a role in real 1-dimensional 
semiconductors [12]. 

Massive Thirring model in 1+1 dimensions can be treated through QISM for both 
bosonic and fermionic field operators [6]. In this article, we shall focus our attention to 
bosonic Massive Thirring (BMT) model. The classical version of BMT model is described 
by the Hamiltonian 

- (0102 + 0200 - 4£0*02020lJ (l.l) 

with the equal time Poisson bracket (PB) relations 

{Mx),Mv)} = {#(*), #(y)} = o, {M^AKv)} = -iS(x-y), 

{0 2 (x), My)} = {$(*), <Pl(y)} = o, {&(*), <P*M} = -i6(x - y) . (1.2) 

It is well known that this BMT model is intimately connected with the derivative nonlinear 
Schrodinger (DNLS) model. In fact, one can generate the Lax operator of BMT model 



by 'fusing' two Lax operators of DNLS model with different spectral parameters [13]. 
The integrability of the classical DNLS model, possesing ultralocal PB structure, can be 
established from the fact that the corresponding monodromy matrix satisfies the classical 
Yang Baxter equation [14]. The quantised version of this DNLS model also preserves the 
integrability property. By applying QISM, the quantum integrability of DNLS model is 
established and the Bethe eigenstates for all conserved quantities have been constructed 
[14,15]. 

In an earlier work by Kulish and Sklyanin [6], the Lax operator and the corresponding 
i?-matrix for the quantum BMT model has been given, though the detailed calculations 
are not being explicitly shown. Moreover, the quantum Yang-Baxter equation (QYBE) 
at the infinite interval limit and hence the corresponding commutation relation between 
the creation and annihilation operators have not been studied. However, it is evident that 
taking the infinite interval limit of the monodromy matrix and corresponding QYBE is 
necessary to get the spectrum for the quantum version of the Hamiltonian (1.1). In this 
context it may be mentioned that, by applying a variant of the QISM [3] which is directly 
applicable to field theoretical models, the quantum DNLS model has been shown to be 
integrable [15,16]. The infinite interval limit of the corresponding QYBE enabled us to 
obtain the spectrum of all the conserved quantities including the Hamiltonian and also 
the two-particle ^-matrix. Therefore, it is interesting to explore the integrability property 
of the quantum BMT model by using the same variant of QISM that we applied for DNLS 
model. In this article our aim is to establish such integrability property of quantum BMT 
model and to obtain the spectrum of all conserved quantities including the Hamiltonian. 

The arrangement of this article is as follows. In Section 2, we consider the classical 
BMT model and evaluate the PB relations among the various elements of the corre- 
sponding monodromy matrix at the infinite interval limit. Using these PB relations, the 
integrability of the classical BMT model can be established in the Liouville sense. In 
this section we also derive the expressions for the classical conserved quantities of BMT 
model. In Section 3, we construct the quantum monodromy matrix of BMT model on a 
finite interval and derive the corresponding QYBE. In Section 4, we consider the infinite 
interval limit of QYBE and obtain the commutation relations among the various elements 



of the corresponding quantum monodromy matrix. Such commutation relations allow us 
to construct exact eigenstates for the quantum conserved quantities of BMT model by 
using the prescription of algebraic Bethe ansatz. In particular we are able to obtain the 
spectrum for the quantum version of the Hamiltonian (1.1). Furthermore we obtain the 
commutation relation between creation and annihilation operators of quasi-particles as- 
sociated with BMT model and find out the S'-matrix of two-body scattering among such 
quasi-particles. In this section we also calculate the binding energy for a iV-soliton state 
of the quantum BMT model. Section 5 is the concluding section. 



2 Integrability of the classical Massive Thirring model 



The classical version of BMT model is described by the Lax operator [6] 

(i{ Pl {x) - p 2 (x)} - ±{A 2 - ^} e{A0t(x) - I0*(:r)} \ 

U(x,\) =i 

V \fa(x) - £&(*) -t{Pl(x) - P 2 (*)} + \{\ 2 - £} ) 

(2.1) 

where p\{x) = 4>l(x)4>i(x), P2{%) = 4>2 ( x ) 4>2(x), A is the spectral parameter and £ is the 
coupling constant of the theory. The bosonic fields <f>i(x), 02 (x) satisfy the PB relations 
(1.2) and vanish at \x\ — > oo limit. The monodromy matrix on finite and infinite intervals 
are defined as 



and 



T* 2 (A) = Vexp r U(x, X)dx (2.2) 

J X\ 

T(X) = lim e(-x 2 ,A)(pexp H U(x,X)dx\e(x 1 ,X) (2.3) 

respectively, where V denotes the path ordering and e(x, A) = e _3 i -j^r^-ix^ 

First, we want to investigate the symmetry properties of the monodromy matrix (2.3). 
It is easy to check that, the Lax operator (2.1) satisfies the relations 

U(x, X)* = KU(x, X*)K, U(x, -A) = K'U(x, X)K' , (2.4a, b) 



where K = \ . , , — 7 _ ) and K ' — I _ . ) . By using these relations, we find 
that the symmetries of the monodromy matrix T(A) (2.3) are given by 

T(X)* = KT{X*)K, T(-X) = K'T(\)K' . (2.5a, 6) 

Due to the relation (2.5a), T(A) can be expressed in a form 

where A is taken as a real parameter. Moreover, by using the symmetry relation (2.5b), it 
is easy to see that o(— A) = a(A) and 6(— A) = —6(A). Therefore, it is sufficient to derive 
the PB relations among the elements of T(X) only for A > 0. 

Next, our aim is to calculate the classical conserved quantities of BMT model by using 
the approach described in Ref 2. From (2.2), one obtains the differential equation followed 
by the monodromy matrix T^ 2 (A) as 

-^T%(\)=U(x 2 ,\)T%(\). (2.7) 

Now, let us decompose the monodromy matrix in the form 

T; i 2 (A) = (l + ^(x 2 ,A))expZ(a; 2 ,x 1 ,A)(l + ^(a; 1 ,A)), (2.8) 

where Z(x2,xi,X) is a diagonal matrix and W(x,X) is a nondiagonal one. The Lax 
operator of the classical BMT model can be expressed as U(x, A) = Ud{x, A) + U n d{x, A), 
where Ud(x,X) is the diagonal part and U n d{x,X) is the non-diagonal part of U(x,X). 
Using the above expression of the Lax operator U(x,X) (2.1), the differential equation 
(2.7) can be decomposed into 

^T = U d + U nd W , ^ - 2U d W - U nd + WU nd W = . (2.9a, 6) 

The structure of the Lax operator (2.1) ensures that W(x 2 , X) and Z(x 2 ,Xi,X) can be 
written in the form 

W(x 2 , A) = -£,w*(x 2 , A)<7+ + w(x 2 , A)<7_ , 
Z(x 2 , x l , A) = z(x 2 , xi, A)<7 3 . 



Substituting eqns (2.6) and (2.8) in the expression (2.3), and using W(x,X) — »• at 
\x\ — > oo limit, one obtains, 

i\ 2 
lna(A) = lim {z(x 2 ,x 1 , X) + — — (x 2 — x^}. 

X2 — >-t-oo zj. 

Xl^ — OO 

Substituting the explicit form of z(x 2 , x±, A) (as obtained by integrating eqn.(2.9a)) to the 
above expression, we get the following form of lna(A): 

/ + OC f + OO nC f + OO 

{0i0i - (f) 2 (f) 2 }dx + it,X / (plwdx-Y <j)* 2 wdx. (2.10) 

-oo J —oo A J —oo 

Next, we expand w(x, A) in inverse powers of A as 

oo 

w(x,\) = Y^ 



°° Wi 



X 2 i +1 

j=0 A 

Using the differential eqn.(2.9b) followed by W(x,X), the expansion coefficient WjS can 
be obtained explicitly in a recursive way. The first few nonzero WjS are given by 

w = -20! ; w 2 = Ai(f> lx + 8^0i (0^0 2 ) + 20 2 . 

Substituting WjS in the expression of lna(A) (2.10), one gets 

lna(A) = ]T-J, 

n=0 A 

where C n s represent an infinite set of conserved quantities. The first two of them are 
explicitly given by 

/+oo 
{0t0i + 0;0 2 }cfe, (2.11a) 

-oo 

/+oo r+oo r+oo 

0*0i x cte + 2£ / {0t0 2 + 0;0i}rfa: + 8e 2 / (PMtyMdx . (2.116) 
-oo J —oo J —oo 

Next we expand w(x, A) in powers of A as 

oo 

w(x,X) = J2tij>? J+1 - 

3=0 

In a similar way as above, using (2.9b), the first few nonzero WjS can be obtained as 
w = -20 2 w 2 = -4i(f) 2x + 8£(0t0l)02 + 20! . 



Correspondingly, eqn.(2.10) yields 

oo 

lna(A) = ^^„A 2n , 

ra=0 

where C n s represent another infinite set of conserved quantities. The first two of them 
are explicitly given by 

r+oo 

C = U { 0101 + 0^02 }cfe , (2.12a) 

J — oo 
„ r+oo r+oo r+oo 

d = Aii / <p* 2 <p 2x dx - 2£ / {0*0 2 + 0*.0 1 }cfe - 8£ 2 / (0*0!)(0;0 2 )rfx. (2.126) 



Now by combining these two sets of conserved quantities, the mass, momentum and 
the Hamiltonian of classical BMT model can be expressed in the following way: 

1 ~ r+oo 

N = --(C -C )= J (0101 + 02> 2 )dx, 



oo 



| „ r+oo 

4£ J —oo 

H = -he, -C 1 )= r [~ i(010 lcc - 0^0 2 J - {0I0 2 + 0^0!} - 4^0^020!] dx. 
4<J J— oo 

Next, we want to derive the PB relations among the elements of T(A) (2.6). We apply 
the equal time PB relations (1.2) between the basic field variables to evaluate the PB 
relations among the elements of the Lax operator (2.1) and find that 

{U(x, \)®U(y, //)} = [r(A, /i), U(x, A) (8) 1 + 1 (8) U(y, /i)] S(x - y) , (2.13) 

where 



r(A, ix) = -£{ t c a 3 <g> a 3 + s c (a + <g> a. + ct_ ® a + ) } (2.14) 

with t c = 2 (\2~- 2) i sC = A 2^ M 2 - Now, by using the eqns.(2.13) and (2.3), one obtains 

{T(A)?T(^)} = r + (A, ^)T(A) ® 7fy) - T(A) ® T(/i)r_(A, //) , (2.15) 

where 

r± = -f (Vcr 3 ® a 3 + s+0+ <8> 0- + sto_ <8> 0+) , 



with s± = ±2i7rA 2 <5(A 2 — fi 2 ). By substituting the symmetric form of T(A) (2.6) to 
eqn.(2.15) and comparing the individual elements in both sides, we obtain 

{a(A),a(/i)} = 0, {a(A),a f (X)} = 0, {6(A), &(//)} = , (2.16a, 6, c) 

{a(A), &(//)} = £ (^rz4) a ( A )%) " 2 ^ A2 ^ 2 " A* 2 ) & (A)a(/i) , (2.16d) 

{a(A), b*(n)} = -$ f ^41 a(A)o» + 2m£A 2 o(A 2 - /i 2 ) 6*(A)a(^) , (2.16e) 



\A 2 -Ai 2 / 
{6(A), &*(//)} = -4mtA 2 o"(A 2 - /i 2 ) |a(A)p 



(2.16/) 



From eqn.(2.16a) it follows that all expansion coefficients occuring in the expansions 
of lna(A) will have vanishing PB relations among themselves. Hence, the following ex- 
pressions will hold true 

for all values of m and n. Since the mass, momentum and the Hamiltonian of the classical 
BMT model has been expressed in terms of the expansion coefficients C n and C n s, all of 
them will have vanishing PB relations among themselves. Thus the integrability property 
of the classical BMT model, described by the Hamiltonian (1.1), is established in the 
Liouville sense. 



3 Commutation relations for the quantum monodromy 
matrix on a finite interval 

By using a version of QISM which is directly applicable to field models [3], in this section 
we shall show that the quantum monodromy matrix of BMT model on a finite interval 
satisfies QYBE. The basic field operators of the quantum BMT model satisfy the following 
equal time commutation relations: 



h(x),My)} = [4>i(x),(j>\(y) 

h(x),fo(y)] = \(f>l(x),(f>l(y) 



0; 
= 0: 



M x )Ai(y)\ = hS(x-y), 

</>2(x),4>l(y) =h6(x-y) 



(3.1) 



and the vacuum state is denned through the relations <pi(x)\0) = 4>2(x)\0) = 0. 

In analogy with the classical Lax operator (2.1), we assume that the quantum Lax 
operator of BMT model is given by 

'fiPi(x) - f 2 p 2 (x) - £ + ^ £\4(x) - i4(x) \ 



U q (x, A) = i 

"' * ' "'' a^ 2 ^ -yivi\*) -ry-2h>2\^) -r -4- — 4^2 



A^i^) - ^02 (x) -^ipi(x) + # 2 p 2 (X) + t ^ 4^ y 



(3.2) 



where pi(x) = 4>\{x)4>i{x) , p 2 (x) = <fi 2 (x)4> 2 (x) and fi,f 2 ,gi, 92 ar e four parameters which 
will be determined later in this section through QYBE. Using the Lax operator (3.2), the 
quantum monodromy matrix on a finite interval is defined as 

Z£ 2 (A) = : Vexp r U q (x, \)dx : , (3.3) 

J X\ 

where the symbol :: denotes the normal ordering of operators. This quantum monodromy 
matrix (3.3) satisfies a differential equation given by 

^-V^X) =:U q (x 2 ,\)T£(\): 

= "^{A 2 - ^sT^(X) +it\<f ) {(x 2 )a + T*>(\) - |4(x 2 )a + T-(A) 
+ i\o-T x ^(\)H*2) - ^_T^(A)0 2 (* 2 ) 
+ tMKx^enT^WM^) - if24(x 2 )e 11 T x x *(\)(j )2 (x 2 ) 
-i 9l ct>\{x 2 )e 22 r^{X)^{x 2 ) +igA{x 2 )e 22 T^{\)^ 2 {x 2 ) , (3.4) 

where en = |(1 + <t 3 ) and e 22 = |(1 — <7 3 ). Now, to apply QISM, we have to find out 
the differential equation satisfied by the product 7^ 2 (A) <g> T^{p). By using the basic 
commutation relations (3.1) and the method of 'extension' [3], we find that the product 
of two monodromy matrices satisfies the following differential equation (detail calculations 
are given in Appendix A): 

d 



dx 2 



(7£(\)®7£(lij) = \£(x 2 ',\,vL)7£(\)®7%(ri\ , (3.5) 



where 



C(x; A, /i) = U q (x, A) <g> 1 + 1 <g> U q (x, p) + C A (x; A, p) , (3.6) 



with 



£ A (x; A,/i) 



/ -hf?Pi(x) 
-hf$p 2 (x) 


-Knh<A{x) 





o \ 





hgifiPi(x) 

+ ^92f2p2(x) 








-h\f i<pi(x) 


-**{>* + & 


hgifipi(x) 


Hi*9i4>i(x) 


-¥Mx) 




+hg2f2p2(x) 


+*f4(x) 




V 


hXgifafa) 





-hgjpi(x) 
-h,g\p 2 {x) 1 



In the expression (3.5), the sign of normal arrangement of operator factors is taken as ::. 
The sign :: , applied to the product of several operator factors (including 0!,02,0i and 
02) 5 ensures the arrangement of all 0{, 2 on the left, and all 0i, 02 on the right, without 
altering the order of the remaining factors. For example, 



;x0 1 2 t 1 0^r; = 0i0 2 xr0 



1^2 ■ 



where X and Y may in general be taken as some functions of the basic field operators. 
Now one can easily check that C(x; A, p) (3.6) follows an equation given by 



R(\, p)C{x\ A, p) = C(x; p, X)R(X, p) , 



(3.7) 



where R(X, p) is a (4 x 4) matrix of the form 

R(X,p) 



(1 0\ 

s(X,p) t(X,p) 

t{X,p) s(X,p) 

\0 1/ 



(3.1 



\ 2 -u? 



with t(A,/i) — A2 A _~^ _ ± , s(A,/i) = |%~^ 2-1 an d q = e~ ta . The above equation (3.7) 
enables us to determine the exact expressions of the parameters f\, f 2 , gi, g 2 , ct in terms 
of the coupling constant £. We obtain : 



K 



sin a, fi = g 2 



& 



-ia/2 



9i = h 



it 



ia/2 



(3.9a, b, c) 



cos a/2' ' " cos a/2 

Using eqns.(3.5) and (3.7), we find that the monodromy matrix (3.3) satisfies QYBE 
given by 

R(X,p)T£(X) ® T£{p) = T£{p) ® T£(X)R(X,p) . (3.10) 



10 



Using the above QYBE (3.10), the commutation relations among all elements of the 
quantum monodromy matrix (3.3) can be obtained easily. 

Eqns. (3.9a,b,c), describing the relations between / 1; / 2 , <?i, # 2 , ol and the coupling 
constant £, provide the necessary conditions for the Lax operator (3.2) to satisfy QYBE 
(3.10). From eqn.(3.9a) we can conclude that, the above method of deriving QYBE for 
quantum BMT model is applicable only when the coupling constant £ lies within the range 
|£| < \- The parameter a has a one-to-one correspondence with the coupling constant 
£ for — | < ct < |. For the purpose of investigating the classical limit of the quantum 
Lax operator (3.2), we take the a — > limit which is equivalent to the % — > limit for a 
fixed value of £. From eqns. (3. 9b, c), it follows that at this limit /i, / 2 — > £ and #i, # 2 — >■ £• 
Hence we find that the quantum Lax operator (3.2) correctly reproduces the classical Lax 
operator (2.1) at h — ► limit. 



4 Algebraic Bethe ansatz for the quantum monodromy 
matrix on an infinite interval 

The quantum monodromy matrix in an infinite interval is defined as 

T(A)= lim e(-x 2 ,X)T^(X)e(x 1 ,X), (4.1) 

X2 — * + oo 
£1^ — OO 

where 7^ 2 (A) is given by eqn.(3.3). Just as in the classical case, the quantum Lax operator 
(3.2) also satisfies the symmetry relations 

U q (x,X)* = KU q (x,X*)K, U q (x,-X) = K , U q (x,X)K', (4.2a, b) 

where K and K' matrices have appeared earlier in eqn.(2.4). Using eqn.(4.2a), the 
quantum monodromy matrix (4.1) can be expressed in a symmetric form given by 

T(A) "U(A) At(A) ) ' (43) 

where A is a real parameter. From eqn.(4.2b), it follows that A(—X) = A(X) and B(—X) = 
—B(X). So it is sufficient to obtain the commutation relations among the elements of the 
quantum monodromy matrix (4.3) only for A > 0. 

11 



Now we aim to obtain the infinite interval limit of the QYBE satisfied by T(X) (4.3). 
To this end, we split the C(x; A,//) matrix (3.6) into two parts: 



C(x; A, /i) = £ (A, //) + d{x\ A, /x) , 



where C (X, ji) is given by 



(--M 2 + ^ 2 ) 



£o(A,/i)= lim £(x;X, /x) 

|a;|— >oo 



+Ka* + 7?) 



-i(A 2 -^ 2 



+ 



i / 1 1 



4VA 2 



f 



-^A/x-fJ i(A 2 -^ 2 



it 1 1 















» f"V2 



V 



IVT2 + 772) / 



4 v A 2 ' /^ 

and C\(x; A, /x) is the field dependent part of C(x; A, /i), which vanishes at x — *■ ±00. From 
eqn.(3.7) we get 

i?(A, fx)e(x; A, /i) = e(x; /x, X)R(X, xx) , (4.4) 

where e(x;X,fi) = e Co<yX '^ x . By using the above mentioned splitting of £(x;A,/x), we 
derive the integral form of differential equation (3.5) as 



re 2 

,T1 



7rf(A) ® ^(/i) = e(x 2 - xi; A, /x) + / dxe(x 2 - x; A, /x) -C^x, A, n)T^(X) <g> Z£(/z): 



From this integral relation it is clear that at the asymptotic limit x±, X2 — ► ±00, Z? 2 (A) <g> 
%c[ 2 {^) ~ * £ { x 2 — X\\ A,/x), which is an oscillatory term. To get rid of this problem, we 
define an operator like 



W{\n)= lim £(-x 2 ;A,/j)T* 2 (A)®T* 2 (/j)£(a;i;A,/i). 

X'2 — > + oo 
X\—> — OO 



(4.5) 



In the above defined operator, the oscillatory nature of T^ 2 (X) ®7^ 2 (/x) has been removed 
and W(X,fi) is perfectly well behaved at the limit X\,Xi — > ±00. By using (3.10) and 
(4.4), it is easy to verify that the operator W(X, /i) (4.5) satisfies an equation given by 



R(X,n)W(X,n) = W(n,X)R(X,ii) 



(4.6) 



12 



The above equation represents the QYBE of BMT model at an infinite interval limit. 

Next, we want to express the QYBE (4.6) directly in terms of the monodromy matrices 
(4.1) defined in an infinite interval. For this purpose, W(X,fi) (4.5) can be rewritten as 



W(X, /i) = C+(A, /x)T(A) <g> T(/i)C_(A, //) 



(4.7) 



where 



C + (A,/i)= lim e(— x; X,fi)E(x; A,/i), C_(A,/i) = lim E(—x;X,fj,)e(x;X,fi), (4.8a, b) 



x— >— oo 



with E(x;\,fi) = e(x, A) ®e(x,n). Substituting the explicit forms of E(x;X,fi) and 
e(x; A, /i) to (4.8a,b), and taking the limits in the principal value sense: linx^ioo P(^-r-) = 
±i7i5(k), we obtain 



C + (\,n) 



/l 0\ 




/I 


0\ 


1 

p+(A,«) 1 

\0 1/ 


, C_(A,/i) = 





Vo 


1 
p-(\,ti) 1 o 
1/ 



(4.9) 



where 



p±(A,/i) 



=F 



=F 



2ihC (A/x + ^ 



1 1 



2 -^T^ + 2 ^( V + ^)' (A2 "" 2 "^ + ^ 



A 2 -p 

2 ^{A/i + ^} 
A 2 - /i 2 - ^ + js =F *e ' 



(4.10) 



Substituting the expression of W(X,/j) (4.7) in eqn.(4.6), we can express this QYBE for 
the infinite interval in the form 



R(X, /i)C+(A, /i)T(A) <g> T(«)C_(A, /i) = C+(//, A)T(«) <g> T(A)C_(/i, A)i?(A, /i) . (4.11) 

By inserting the explicit forms of i?(A,/x) (3.8), C±(X, fi) (4.9), and T{X) (4.3) to the 
above QYBE (4.11) and comparing the matrix elements from both sides of this equation, 
we obtain the following commutation relations: 



[A(X),A(f,)] = 0, \A(X),A\p)]=0, [B(\),B(»)]=0 



A(X)B\n) 



fi 2 q — X 2 q~ 1 
« 2 — A 2 — ie 
H 2 q — X 2 q~ l 



B\n)A{X) 



/i 2 -A 



., B\y)A{\) - 2nTiiXp8(X 2 - p?)B\X)A{p) 



(4.12a, 6, c) 



(4.12a 7 ) 



13 



Li /\ 6c 

= ^KKl W )£fc) - 27r^A/!<5(A 2 _ ^)A{n)B{\) , (4.12e) 

B(n)B\\) = r(A, fi)B\X)B{fi) + 47r^A/i5(A 2 - fi 2 )A\X)A{\) , (4.12/) 

where 

8^ 2 AV 4^ 2 {A/i+^} 2 



r(A,//) 



1 



(A 2 .^ (A 2 _^2 _ ^ + J^ _, e)(A 2 _^2 _ £ + ±_ + fc) 

It is interesting to note that, for the case A ^ /J, eqn.(4.12f) gives [-B(A), -B^^u)] 7^ 0, 
whereas from eqn.(2.16f), one obtains that {6(A), &*(//)} = for A^/i. Thus the corre- 
spondence principle is not manifest here in a straightforward manner. However the h — ► 
limit of r(A, /i), gives the correct classical counterpart of the commutation relation (4.12f ). 
Due to eqn.(4.12a), all the operator valued coefficients occuring in the expansion of 
lnA(A) will commute among themselves. As a consequence the BMT model described 
by the Lax operator (3.2) turn out to be a quantum integrable system. By applying 
the method of algebraic Bethe ansatz, one can also construct the exact eigenstates for 
all commuting operators which are generated through the expansion of lnA(A). With 
the help of eqn.(4.1), it is easy to find that A(A)|0) = |0). By using this relation and 
eqn.(4.12d), it can be shown that 

A(X) \fx u /x 2 , • • • , fx N ) = ]j{ r 2 _ , |/ii, /x 2 , • • • , Hn) , (4.13) 

where fj,jS are all distinct real or complex numbers and |/xi, ^2, • • • , /j>n) = B^(ni)B^(n 2 ) ■ ■ ■ 
B '(jUjv)|0) represents a Bethe eigenstate. Using the commutation relation (4.12f) one 
can also calculate the norm of the eigenstates B^(fii)B^(^ 2 ) • • • B t(//jv)|0). However, the 
commutation relation (4.12f) contains product of singular functions (A 2 — fi 2 — -X + \ — 
ie) _1 (A 2 — fi 2 — -h + \ + ie)^ 1 , which is undefined at the limit A — ► /i. As a result, 
eigenstates like B*(hi)B\ij,2) • • 'B'(ij,n)\0) are not normalised on the ^-function. To 
solve this problem, we consider a reflection operator given by 

i?t(A) = Bt(A)(At(A))' 1 (4.14) 
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and its adjoint R(X). By using eqns.(4.12a-f), we find that such reflection operators satisfy 
well defined commutation relations like 

R\\)R^fi) = S-\\, /i) rf(fi)R j (\) , 

R(X)R(fi) = S-\X, /i) R(/j)R(\) , 

R{X)R\fi) = S{\, /i) R\fi)R{\) + AnhX 2 5(X 2 - /1 2 ) , (4.15) 

where 

g(A,/0= xXT ^C ■ ( 4 - 16 ) 

X 2 q L — fi 2 q 
The iS(A, /i) defined above represents the nontrivial S'-matrix element of two-body scat- 
tering among the related quasi-particles. We find that this S(X, //) satisfies the following 
conditions: 

S-\\, f i) = S( t i,\) = S*(\,fi), (4.17) 

and remains nonsingular at the limit A — ► \x. Consequently, the action of the operators 
like R'(X) on the vacuum would produce well defined states which can be normalised on 
the ^-function. 

The point to be noted here is that in eqn.(4.13), the eigenvalues of A(X) are in general 
complex. To get real eigenvalues, we define a new operator In A(X) through the relation 
lnA(A) = In A(Ae~2~) and expand this operator in inverse powers of A: 

MA) = E§- (4.18) 

n=0 A 

Using eqns.(4.13) and (4.18), it is easy to see that C n s satisfy eigenvalue equations like 

C n |a*1, A*2, • • • , A*jv) , = Xn |/^1, A*2, ' ' ' , Hn) , 

where the first few XnS are explicitly given by 

N N 

Xo = aN, xi = 2sina^/i 2 , X2 = sin 2a J2 1$ ■ ( 4 - 19 ) 

3=1 3=1 

It may be noted that these eigenvalues are all real when [XjS are taken as real numbers. 
Next we expand the operator lnA(A) in powers of A as 

oo 

lni(A) = ^C„A 2n , (4.20) 

n=0 
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and by using (4.13) we obtain 

C n |/ii, /i2, • • • , /ijv) , — Xn \fJ>l, A*2, ■ • • i A*jv) • 

The first few x„s are explicitly given by 

N i AT i 

^ = - a iV, xi = -2sina^ — , £2 = - sin 2a £) — . (4.21) 

In analogy with the classical case, one can now define the momentum and Hamiltonian 
of the quantum BMT model as 

P = ~(C 1 + C 1 ) 1 H = -±(C 1 -C 1 ). 

By using (4.19) and (4.21), the eigenvalue equations corresponding to the above momen- 
tum and Hamiltonian are obtained as 



1 N ( 1 \ 

P|/!l, /i 2 , • • • , //jv) = - J2 [$ 2 

2 i=l \ l*jJ 



1/^1)^2, • • • , fJ'N) 



1 N ( 1 \ 

W|/xi, /x 2 , • • • , /Mr) = - J2 [$ + T I lA*i» A*2, • • • , /ijv) • (4-22) 

In the above expressions, /i^s are taken as real numbers and |//i,// 2 , • • • , (J>n) represents 
a scattering state. Now to construct quantum iV-soliton states of BMT model, complex 
values of fij can be chosen in such a way so that the eigenvalues corresponding to different 
expansion coefficients of In /1(A) still remains real. Such a choice is given by 



Hj = \x exp 



. {N+l . 

V 2 J 



(4.23) 



where jj, is a real parameter and j G [1, 2, • • • , N]. For the above choice of /ij, eqn.(4.13) 
takes the form 

_\ 2 _ 2 N+l 

A(A)|/zi,/i 2 ,---,^v) = g~ JV ( A2 _ 2 _j V+ i )l^i>/ i 2,---^Jv)- (4.24) 

Consequently, the energy eigenvalue equation corresponding to the quantum iV-soliton 
state can be obtained as Tt\/J>i, /i 2 , • • • , A* at) = E\/j,i, /i 2 , • • • , jxn) , where 

1 / o 1 \ sinaA^ 



2 \ ji z I ana 
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Thus we find that quantum N-soliton states can be constructed for BMT model for N > 1. 

Now we assume a particular value of the coupling constant £ given by £ = — sin a = 

— sin(^^), where m and n are nonzero integers which do not have any common factor. 

From eqn.(4.23), we obtain \x^ = fij+n for the above choice of £. Since all the ^s have to 

be distinct, we get N < n as a restriction on the number of quasi-particles that can form 

a quantum soliton state for BMT model when £ = — sin( 2iL! ^). 

Next we aim to calculate the binding energy for a iV-soliton state of quantum BMT 

model. Substituting the expression of fij (4.23) to the first relation in eqn.(4.22), the 

momentum eigenvalue of a iV-soliton state is obtained as 

„ l,o 1 . sin aN 
2 fi z sin a 

It is intersting to observe that the energy (4.25) and the momentum eigenvalue (4.26) of 

a iV-soliton state satisfy the dispersion relation E 2 = P 2 + m 2 , where m = sin aN . To 

calculate binding energy we assume that the momentum P (4.26) is equally distributed 

among N number of single-particle scattering states. The real wave number associated 

with each of these single particle states is denoted by /i - With the help of eqns.(4.22) 

and (4.26), we find that 

9 1 / 9 1 N sin aN , , . 

IM --= /i 2 -- ) w —■ 4.27 

Using eqn.(4.22), the total energy for iV number of such single particle states is obtained 

as 

1 
r,, N ( 2 !\ N \ ( 2 1 V sin 2 aN 1 5 lt , 

E= j{^ + J4) = j{^-7)^^ + T (428) 

Subtracting E (4.25) from E' (4.28), we obtain the binding energy of the quantum 
iV-soliton state as 

E B (a,N) = E'-E 

sin^ aN I 1 / 9 1 \ sin aiV 

- ) Ar2 • 2 + 4 ? " o U + -2 — • 4 ' 29 

i 2 J N 2 sin 2 a J 2 \ /i 2 / sina 

Note that the above expression of Es{a,N) remains invariant under the transformation 
a — ► —a. So it is sufficient to analyse the nature of binding energy within the range 
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< a < f . Now, for E B (a,N) to represent the energy of a real bound state, E' has to 
be greater than E. Since E' (4.28) is always positive, it is evident that E' > E for E < 0. 
So we will restrict our attention only for the case E > 0, when the condition E' > E is 
equivalent to -E' > E 2 . Substituting the explicit expressions for E' (4.28) and E (4.25), 
the above condition takes the form 

N sin a > sin aN . (4.30) 

Substituting N = 2 in (4.30), we get the trivial inequality 1 > cos a for a > 0. So the 
condition (4.30) is satisfied for N = 2 case within our chosen range of a. By using the 
method of induction, we can easily prove that the condition (4.30) is valid for arbitrary 
values of N. Thus we get an iV-soliton bound state when a lies in the range < \a\ < |. 

5 Concluding Remarks 

In this article we consider the classical Lax opeartor of BMT model and obtain the PB 
relations among various elements of the classical monodromy matrix at the infinite interval 
limit. By using these PB relations, the classical integrability of BMT model is established 
in the Liouville sense. We also calculate the classical conserved quantities of BMT model. 
Next, we quantise the Lax opeartor of BMT model. By using a variant of QISM, that can 
be directly applied to the field theoretic models, we obtain the QYBE for the quantum 
monodromy matrix at a finite interval. This QYBE enables us to determine the various 
parameters of the quantum Lax operator in terms of the coupling constant £. Then we 
take the infinite interval limit of this QYBE and derive all possible commutation relations 
among the various elements of the corresponding quantum monodromy matrix. These 
commutation relations enable us to establish the quantum integrability of BMT model 
and also to construct the exact eigenstates for the quantum version of the Hamiltonian 
(1.1) as well as other conserved quantities by using algebraic Bethe ansatz. We also 
obtain the commutation relation between creation and annihilation operators associated 
with quasi-particles of BMT model and find out the ^-matrix for two body scattering. 

In this context, we consider the BMT model with some special values of coupling 
constant given by £ = — sin a = — sin(^p), where m and n are nonzero integers with no 
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common factor. It turns out that the number of quasi-particles, which form a bound state 
for such quantum BMT model, cannot exceed the value of n. We have also derived the 
exact expression of binding energy for a iV-soliton state of quantum BMT model. The 
binding energy turns out to be positive for all allowed values of a. 

The commutation relation between creation and annihilation operators will play an 
important role in the future study, since by using it one might be able to calculate the 
norm of Bethe eigenstates and various correlation functions of the BMT model. In future, 
we would also like to obtain the quantum conserved quantities of BMT model in terms 
of the field operators by using a method which was used earlier in the case of nonlinear 
Schrodinger model [17] and DNLS model [16]. 
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Appendix A 

Here we give the details of deriving eqn.(3.5). Direct attempt to calculate 



a 

3x2 



(Z^ 2 (A) ® Xnif 1 )) by using eqn.(3.4), leads to indeterminate expressions of the form 



T£(\),4{x 2 ) 



and 



7? 2 (A), 2 (x 2 ) • To avoid this problem by using the method of 
extension [3], we shift the upper limit of the monodromy matrix 7^ 2 (A) by a small amount 
e and take e — > limit only after differentiating the product T x X2+e (X)®)T X2 (n) with respect 
to x 2 . So, using eqn.(3.4), we obtain 
d 



dx 2 
where 



T x ^(\)®T**(n) 



I ( U q {x 2 + e; A) <g> 1+ 1 <g> U q {x 2 ; fi) ) T x X2+e (\) ® T x X2 (n) i 
+ K + + if_, (Al) 



^ 



K_ 



(A)® 2 (x 2 + e),T a f» 



= ^ T*^(\),4(x 2 ) w+TZM-S T x ^(X),4(x 2 ) ®a + T^ 2 (/i) 

L J yU L J 

+ i/i [T£ + *(X),4(x 2 )] ^en^^ffe) 

- i/ 2 [^ 2+e (A),^(x 2 )] ® eil 7^(/x)^(x 2 ) 

- x^i [7^ 2+e (A)>i(x 2 )] QevTZMfa (x 2 ) 
+ ig 2 [T x X2+ *(\), 2 (x 2 )] ® e 22 T£ {n)ct> 2 {x 2 ) , 
= iAa_T- +e (A) ® [0!(x 2 + 6),T-(/i)] - ia_7^ 
+ i/iri(x 2 + c)en7^ +e (A)® [^(x. + e),^ 2 ^) 

- */ 2 2 (x 2 + e) eil T^ 2+e (A) ® [0 2 (x 2 + e), 7£» 

- ^i0l(^2 + e)e 22 T x x ^(X) ® [^(x, + e),T^(^) 
+ ^4(^2 + e^T^A) ® [0 2 (x 2 + e), ^(/x) 

Now we consider the case, e > 0. Since 0i(x 2 + e) and 2 (x 2 + e) commute with 
'i(x), 0J(x), 2 (x), 2 (x) for all x lying within xi and x 2 , we get 0i(x 2 + e),T x 2 (n) = 
</>2(x 2 + e), Z^ 2 (/i) = 0. Thus we can conclude that for a positive e, K_ = . So we have 

to calculate only the nontrivial commutators Z C X2+<E (A), 0i(x 2 } 

appearing in the expression of K + . 



and 



r^ 2 ^(A),0 2 (x 2 



First let us calculate the commutator 



7^ 2+e (A), 0i (x 2 ) . For this purpose, we consider 
a 'transformation' Q, which replaces the classical variables ip(x) and ip*(x) by quantum 
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operators ip(x) and i]y{x') respectively. Next we use a correspondence principle [3], 



T x ^(X),4(x 2 )} =ih:n{TZ +e (q;\),<j>* 1 (x 2 )} :, (A2) 

where T^ +e (q; A) represents a classical monodromy matrix given by 

T^ +e (q; X) = V exp / U q (x, \)dx , 

and U q (x,X) = fl~ 1 U q (x,X). By using the fundamental PB relations (1.2), it is easy to 
find that 

rx 2 +e 

{T-+%A),0J(:r 2 )} = / dxT^{ q -X){U q {x,X)A\{x 2 )}T:^X) 

JX\ 

= Tlf\q- X)(fMx 2 )e u - 9l ^(x 2 )e 22 + Xa.)T^(q; A) . 



Taking e — > limit of the above expression and substituting it in (A. 2), we obtain 



V? +e W,<ilM] = ^(fi4(x 2 )e 11 -g 1 (f ) \(x 2 )e 22 + Xa) T£{X) 



lim 



Next we have to calculate the commutator 



(A3) 



7^ 2+e (A), (j) 2 (x 2 ) • Using the same corre- 



spondence principle as before and finally taking the e — > limit one obtains, 



lim 



1 



^? +e (A), 4(x 2 ) ]=ih( -f2<A(x 2 )en + 924>\{x 2 )e 22 - -a. T%(\) . (A4) 



A 



Taking the e — > limit of eqn.(Al) and using (A3) and (A4), we finally obtain the 
differential equation (3.5). Note that, instead of e > 0, we could have chosen e < in 
eqn.(Al). In that case only the commutators 4>\(x 2 + e), 7^ 2 (/i) and 4> 2 (x 2 + e),T x x ^(/j,) 
give nontrivial contributions. However, by repeating similar steps as outlined above and 
finally taking the e — > limit, we would have obtained the same differential equation 
(3.5). 
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